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ABSTRACT 

We present new results on the shock around the southwest radio lobe of Centaurus A using 
data from the Chandra Very Large Programme observations (740 ks total observing time). 
The X-ray spectrum of the emission around the outer southwestern edge of the lobe is well 
described by a single power-law model with Galactic absorption - thermal models are strongly 
disfavoured, except in the region closest to the nucleus. We conclude that a significant fraction 
of the X-ray emission around the southwest part of the lobe is synchrotron, not thermal. We 
infer that in the region where the shock is strongest and the ambient gas density lowest, the 
inflation of the lobe is accelerating particles to X-ray synchrotron emitting energies, similar 
to supernova remnants such as SN1006. This interpretation resolves a problem of our earlier, 
purely thermal, interpretation for this emission, namely that the density compression across 
the shock was required to be much larger than the theoretically expected factor of 4. We 
describe a self-consistent model for the lobe dynamics and shock properties using the shell of 
thermal emission to the north of the lobe to estimate the lobe pressure. Based on this model, we 
estimate that the lobe is expanding to the southwest with a velocity of ^2600 km s _1 , roughly 
Mach 8 relative to the ambient medium. We discuss the spatial variation of spectral index 
across the shock region, concluding that our observations constrain j max for the accelerated 
particles to be ^10 8 at the strongest part of the shock, consistent with expectations from 
diffusive shock acceleration theory. Finally, we consider the implications of these results for 
the production of ultra-high energy cosmic rays (UHECRs) and Te V emission from Centaurus 
A, concluding that the shock front region is unlikely to be a significant source of UHECRs, 
but that TeV emission from this region is expected at levels comparable to current limits at 
TeV energies, for plausible assumed magnetic field strengths. 

Key words: galaxies: active - galaxies: elliptical and lenticular - galaxies:individual (Cen 
A) - radio continuum: galaxies - X-rays: galaxies - shock waves 



1 INTRODUCTION 
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The mechanisms by which radio-loud active galaxies transfer en- 
ergy to their surroundings are of considerable interest, due to the 
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likely role of such energy input in counteracting gas cooling in the 
central regions of galaxies, galaxy groups and clusters. Shocks are 
particularly important, as they provide a means of increasing the 
entropy of the surrounding gas. There are now a number of exam- 
ples of shocks associated with radio galaxies, which include weak 
shock structures in the intr acluster medium (ICM) of some of the 
neare st galaxy clusters (e.g. lFabian et al.l200 3. 2006; Form anet al.l 
2005) generally on scales larger than the radio lobes, and a small 
number of examples of strong shocks associated with g alaxy-scale 
radio lobes in comparatively isolated environments (e.g. lKraft et al.l 
l200l l2007d : ICroston et alj2007l : Ijetha et alj|2008l) . 

The best-studied example of a strong shock associated with 
a radio galaxy is the prominent shell of bright X-ray emission 
surrounding the south-west inner lobe of Centaurus A (Cen A), 
detected w i th Ch andra and XMM-Newton, as first reported by 
iKraft et alj d2003l) . As Cen A (NGC5128) is the nearest radio 
galaxy (Dl = 3.7 Mpc, see below), it is an ideal laboratory for in- 
vestigating the physics of radio-lobe shocks and their effects on the 
galaxy interstellar medium (ISM). Kraft et al. (2003) interpreted the 
bright X-ray emission as thermal in origin, based on acceptable fits 
to a thermal spectrum, and a lack of radio emission at the location 
of the shock front (as might be expected if the X-ray emission was 
associated with high-energy particle acceleration); however, their 
analysis raised questions about the shock hydrodynamics, as the 
apparent density contrast between shocked and unshocked material 
was significantly higher than the expected va lue of 4 for a stron g 
shock (where the adiabatic index is 5/3). In IKraft et all d2007al) . 
deeper Chandra observations were used to resolve the structure of 
the shock front in more detail, and to investigate the shock structure 
as a function of position. They developed a more complex two-fluid 
model for the shocked material, which accommodated the electrons 
and ions in the shocked material being out of equilibrium. However, 
that work did not address the issue of the density contrast between 
pre-shocked and post-shocked gas. 

Here we present a new analysis of the Cen A shock properties, 
based on very deep Chandra observations of the source as part of a 
Very Large Programme carried out in Chandra Cycle 9. A total of 
~ 740 ks of Chandra data on Cen A is now available by combin- 
ing all of the available Chandra data. This rich dataset has already 
been used to carry out detailed studie s of the X-ray binary popula- 
tion in NGC 5128 djordan et al.l2007l;lsivakoff et al.l2008h. particle 
accel eration in the X-rav jet dHardcastle et alj|2007l ; IWorrall et al\ 
l2008h . and of the galaxy ISM dKraft et alj|2008l) . Our analysis in 
this paper of the new deep observations of the Cen A shell revises 
some of the conclusions of earlier work by Kraft et al. 

In Section [2] we summarize the data reduction and analysis, 
in Section[3]we discuss the results of our new spatial and spectral 
analysis of the X-ray shell, and in Section|4]we compare our results 
with previous work, present a new interpretation for the majority of 
the shell X-ray emission, and discuss the implications of that inter- 
pretation. We adopt a distance of 3.7 Mpc to Cen A (the average of 
five d istance estimates to Cen A - see Section 6 of iFerrarese et al.l 
|2007|) . giving an angular scale of 1 arcsec = 17.9 pc. 
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Figure 1. The shell of X-ray emission surrounding the SW lobe, from the 
merged Chandra data from all observations to date, in the 0.4 - 2.5 keV 
energy range. The image has been lightly smoothed with a Gaussian kernel 
of FWHM 2.3 arcsec. Radio contours at 1.4 GHz are overlaid with contour 
levels 9.0 X (1, 4, 16, 64) mjy beam -1 . The colour scale is linear in units 
of ACIS counts/pix. 

entire shell region is ~ 507 ks. The inner parts of the shell have a 
higher total exposure time, but for this analysis we chose to include 
only the datasets in which the entire shell falls within the detector 
field of view. The datasets were all processed uniformly using CIAO 
version 3.4 and CALDB version 3.3.0. 1 . Furt her d etails of the data 
proces sing are given in lSivakoff et alj d2008h and lHardcastle et all 
d2007l) . Spectral fitting was carried out using XSPEC. 

Since the primary beam of the VLA limits the sensitivity of 
maps of Cen A at high frequencies and large distances from the 
pointing center, we made a new radio map with an effective fre- 
quency of 1.5 GHz from archival VLA datasets. We combined data 
from observations with the VLA in its A configuration (AB257: 
1983 Oct 30), B configuration (AC175: 1986 Aug 11), BnC con- 
figuration (AC175: 1986 Oct 03) and CnD configuration (AB587: 
1991 Feb 16) to give a map that had the highest possible resolu- 
tion while sampling the large-scale structure to the extent possi- 
ble with the VLA at this frequency. Calibration and data reduction 
were carried out with AIPS in the standard way. The resulting map 
has a resolution of 6.1" x 1.9" (FWHM of major and minor axes 
of elliptical Gaussian: the major axis is closely aligned with the 
north-south axis) and is dynamic-range limited by artefacts around 
the bright core, with an off-source r.m.s. around 1.8 mjy beam . 
Images presented in this paper are corrected for attenuation due to 
the primary beam of the VLA. 



2 DATA PREPARATION AND ANALYSIS 

Cen A has now been observed a total of 10 times using Chandra's 
AXAF CCD Imaging Spectrometer (ACIS) instrument. The south- 
west radio lobe and surrounding X-ray shell are entirely within the 
detector area for 7 of these observations (ObsIDs 0962, 2978, 3965, 
7797, 7800, 8489 and 8490). The effective on-source time for the 



3 RESULTS 

Fig. [T] shows a Gaussian-smoothed image of the X-ray emission 
surrounding the south-west radio lobe with radio contours overlaid. 
It can be seen that the shell is sharply edge-brightened. We interpret 
the outermost edge of the X-ray emission as the shock front. This 
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Figure 3. Source and background regions used for the spectral and spatial analysis. Top left: Region 1, which includes the brightest part of the shell emission; 
top right: Regions 2 and 3: the region of X-ray shell emission exterior to the radio lobe (bottom right of the image) and the region of bright emission at the 
north-east edge of the lobe, respectively; bottom left: zoomed-in view of Region 3, the north-east edge of the SW shell; bottom right: regions used for surface 
brightness profile in Fig. [2] Chip gaps and out of field-of-view regions were masked out of background spectra on an observation by observation basis as 
necessary, and point sources were excluded from all spectra. 
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Figure 2. X-ray (solid line) and radio (dashed line) surface brightness pro- 
files across the shock front (the extraction regions are shown in Fig. [3). Ar- 
row symbols indicate positions at which the radio profile is based on (3ct) 
upper limits. The profiles have been renormalised for ease of comparison. 
The dotted line indicates the position of the shock front as determined from 
the X-ray data. The bin width is significantly larger than the size of the radio 
beam, so that all data points are independent. 



lies ~ 16 arcsec from the outer edge of the radio lobe, which is well 
defined in our radio images, as demonstrated by the radio and X-ray 
surface brightness profiles shown in Fig.(2] extracted for the regions 
shown in the bottom right panel of Fig. [3] We interpret the sharp 
radio lobe boundary to be the contact discontinuity. Shell emission 
is seen throughout the lobe region, including a bright region on the 
northern edge of the shell, a d istance of 1 . 4 arcmi n E of the nucleus, 
which was first discussed bv lKraft et alj d2007al) . The fact that the 
X-ray emission clearly extends a significant distance beyond the 
edge of the radio lobe (well defined in our high resolution radio 
maps) and is edge-brightened (see Fig.|2j strongly supports a model 
in which the X-ray emission originates in a shell surrounding the 
radio lobes, rather than originating in a fluid contained within them. 

We initially carried out spectral analysis for several large re- 
gions to determine the global properties of the shell X-ray emission. 
The spectral extraction regions are shown in Fig. [3] Spectra were 
extracted from all seven datasets using the CIAO task specextract, 
which also calculates appropriate response files. The background 
regions used are also shown in Fig. [3] For the background spec- 
tra, additional exclusion regions were included on an observation- 
by-observation basis to mask any out-of-field-of-view regions and 
chip gaps. Point sources were also excluded from all of the ex- 
tracted spectra. The spectra were fitted in the 0.4 - 3.0 keV energy 
range to minimize background contamination from the wings o f 
the PSF from the very bright central AGN (e.g. lEvans et alj|2004[) . 
which dominate above 3.0 keV across the entire field-of-view (e.g. 
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iKraft et al.l2 007a). This choice of energy range also minimizes any 
contamination from readout streaks. In the few cases where read- 
out streaks pass through regions of interest (e.g. ObsID 7800 for 
Region 3, and ObsIDs 7797 and 8490 for a small number of inner 
regions discussed in Section [3~3t we checked that excluding the 
datasets that could be contaminated did not affect the spectral re- 
sults. For each region we carried out joint fits to the spectra from 
all epochs, and we assumed a fixed hydrogen c olumn density at the 
Galac tic value of N H = 8.4 x 10 20 cm -2 jPickev & Lockmanl 
1990), except where otherwise specified. This is reasonable in gen- 
eral, as our spectral extraction regions are mainly well away from 
the central region where varying absorption from gas associated 
with the dust lane needs to be taken into account. For each of 
the regions we initially fitted five s pectral models: an APEC (As- 
trophysical Plasma Emission Code: ISmifh et alj |2001 ) model with 
free abundance (Model I), an APEC model with abundance fixed 
at 0.15^0 (Model II), an APEC model with abundance fixed at 
0.5Zq (Model III), a single power-law model (Model IV) and a 
model consisting of an APEC (Z = 0.5) plus a power-law com- 
ponent (Model V). The abundance values for Models II and III 
were chosen to sample the range of abundances from the low val- 
ues found in free- abundance fits to m ore typical values for galaxy 
atmospheres (e.g. iDavid el a/.ll2006t) . Complete results of all the 
spectral fits are given in Tabled Quoted errors are la for one or 
two interesting parameters except where otherwise specified. 

3.1 The bright outer shell emission 

We first considered a region encompassing all of the bright shell 
emission (Region 1 in Fig. [3j, with point sources excluded. While 
the spectral model parameters vary somewhat across this large re- 
gion, the general features of the model fits described below are 
present in all the fits from smaller regions when Region 1 is subdi- 
vided. For this region we find that Models IV and V (single power 
law and 2-component power-law plus APEC, respectively) give the 
best fits to the spectra. Models I - III (thermal models) give much 
poorer fits. Model I, the free abundance fit, has the lowest \ 2 °f 
all of the thermal models; however, the abundance is constrained 
to be extremely low (a 3<r upper limit on Z of 0.007 Zq), and the 
fitting statistic is still significantly poorer than for the power-law fit 
(Model IV). For Model V (APEC plus power law), the inclusion 
of a thermal component does not significantly improve the fit when 
compared to Model IV, and the thermal model contributes an unab- 
sorbed 0.4 - 3.0 keV flux of ~ 1.6 x 10~ 14 erg cm -2 s -1 , whereas 
the power-law contributes a flux of ~ 4.7 x 10 -13 erg cm -2 s _1 , 
a factor of 30 times higher. The spectrum for Region 1 is shown in 
Fig. [4] with each of the five models overplotted. The fitting statis- 
tics for this region are fairly poor for all five models. This is likely 
to be due to radial variation in the spectral properties of the emis- 
sion, possibly including variation in the contribution from thermal 
emission. We explore this spatial variation in Section [33l 

To summarize, we find that a single power law is a better de- 
scription of the spectrum than any thermal models. In addition, the 
fitted temperature values for thermal models are strongly depen- 
dent on the adopted abundance. T he shallower observations anal- 
ysed bv lKraft et all ( l200i 12007 j) meant that thermal models pro- 
vided an adequate fit to the data, but with our deeper data, we are 
now able to rule out a simple, single-temperature APEC model at 
high confidence. In particular, the thermal fits have large residuals 
at soft energies. Fitting a two- temperature APEC model can miti- 
gate this problem to some extent: we find best-fitting temperatures 
of 0.28 ± 0.02 keV and 3.90 ± 0.26 keV, assuming 0.5 times solar 



abundances. This fit has \ 2 of 1065 for 742 d.o.f, which is a sig- 
nificant improvement on the single temperature models for realistic 
abundances, but still considerably poorer than the power law fit. 

In order to firmly rule out a thermal origin for the shell emis- 
sion, we also considered the limits on the 6.7 keV line from ionized 
Fe, and the Fe L complex, for assumed gas temperatures of 1.5 and 
3 keV (spanning the range seen in the thermal spectral fits). We 
found that the limits on the 6.7 keV line do not provide strong con- 
straints, but the predicted fluxes in the Fe L complex at 1.17 keV 
for an abundance of 0. 15Z0 are 1.6 x 10~ 5 and 4 x 10~ 6 photons 
cm~ 2 s _1 , respectively, whereas the 3<j upper limit from our joint 
fits (using the linewidth and energy of a broad Gaussian fit to the 
complex in a faked spectrum of the same parameters) is ~ 2 x 10~ 6 
photons cm -2 s _1 . Hence there is no evidence for the line emis- 
sion expected from a thermal plasma with realistic temperature and 
abundance values. We cannot rule out a model in which turbulent 
broadening has erased the lines from the spectrum, but this would 
require turbulent velocities at or well above the sound speed of the 
gas. 

We also examined separately the brightest region of shell 
emission at the south-west rim, ~ 5.6 arcmin to the SW of the 
nucleus (Region 2 in Fig. [3}, which is external to the radio lobes 
(see Fig. [2}, to construct a broad-band SED for the shell without 
contamination from any possible radio lobe emission (see below). 
Details of the spectral fits are in TableQ] The overall conclusions of 
the spectral analysis for this region are very similar to those for Re- 
gion 1 ; this is, the spectrum is much better fitted by a non-thermal 
power-law component rather than thermal emission. This is as ex- 
pected given that this region contributes a significant fraction of the 
flux in Region 1 . 

We therefore conclude that the shell emission is likely to be 
non-thermal, in contrast to the conclusions of earlier work based 
on poorer quality spectra. 

3.2 The north-east edge of the shell 

As shown in Fig.[T] there is a small region of brighter shell emission 
at the northern edge of the SW shell close the AGN (~ 1.4 arcmin 
east of the nucleus: Region 3 in Fig.[3). This region w as shown to 
be asso ciated with the shell and examined in detail by IKraft et al.l 
d2007d) , who concluded that the emission was thermal, and signifi- 
cantly cooler than at the outer, SW edge of the shell (0.6 — 0.8 keV 
for their SI and S2 regions for an assumed abundance of 0.5Zq). 
Since the line emission is expected to be stronger at these temper- 
atures, this region can potentially provide useful constraints on the 
shell X-ray emission. We t herefore examine d spectra for Region 3 
(Fig.|3j, which is similar to lKraft et all J2007ah 's S 1 region. Results 
are listed in Table[T] 

In contrast to the situation in Region 1 and sub-regions of the 
outer shell, we find that a thermal model is strongly favoured in 
this region. With a free abundance fit, we obtain a best-fitting tem- 
perature of 0.91 ± 0.04 keV and an abundance of 0.17±^;^| Z 
with x 2 = 137 for 97 d.o.f. The Fe L complex is very prominent in 
the spectrum, which cannot be fitted adequately with a power-law 
model. The inclusion of a power-law component to the model, in 
addition to the APEC, does not improve the fit. Hence we conclude 
that the X-ray emission from this region is unambiguously thermal. 
The best-fitting abundance is significantly lower than that expected 
from the galaxy ISM, which is likely to be due to the assumption 
of a single temperature (e.g. lBuotell200oT) . As our spectral regions 
are likely to include multi-temperature gas, we cannot obtain a re- 
liable estimate of the metallicity from the X-ray data. We therefore 
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Figure 4. Spectral fits for Region 1. Top left: Model I (free-abundance APEC), middle: Model II (APEC with Z = 0.15Zq), right: Model III (APEC with 
Z = 0.5Zq), bottom left: Model IV (power law), right: Model V (power law + fixed-abundance APEC). The data are in red and the fitted models are in 
green. The y-axis scale is different for Model V due to the inclusion of the thermal component, which has a low flux. 



chose to adopt a metallicity of Z — 0.5Zq, more appropriate for a 
galaxy halo. The single temperature fits also showed strong resid- 
uals at ~ 1.35 keV, likely to be due to Mg XI, and so we chose 
to include a narrow Gaussian at this energy. The best fitting tem- 
perature for an APEC model with fixed abundance plus 1.35 keV 
Gaussian is 0.95± o ;^ keV, with \ 2 = 145 for 97 d.o.f. We also 
investigated whether additional Silicon features are present in the 
spectrum as would be expected if the gas is overabundant in alpha 
elements. We found a significant improvement in the fit if a Gaus- 
sian is added at 1.865 keV, corresponding to Si XIII, and a marginal 
improvement including a Gaussian at 2.006 keV, corresponding to 
Si XIV. Using the VAPEC model with all other elements fixed at 
an abundance of 0.5 times solar (which gives a very similar fit), we 
find a magnesium abundance of 2.4^0 4 times solar. If the silicon 
abundance is allowed to vary, an abundance of 1.2 jl^'l times so- 
lar is obtained. The data are therefore consistent with a model in 
which the shocked gas in this region is alpha-enriched. It is also 
possible that the unexpectedly strong magnesium and silicon lines 
are indicative of a non-equilibrium ionization state. 

It is likely that there is complicated temperature structure in 
the shell that we cannot resolve, and as the gas in this region is as- 
sumed to have passed through the shock front recently, it is likely 
that the ionization balance has been affected by the shock. We fitted 
non-equilibrium ionization models to the spectra for Region 3 (the 
nei model in XSPEC), and found a fit of similar quality to that for 
the APEC model with a Gaussian at 1.35 keV, with a very similar 
temperature and normalisation to our adopted model. Our conclu- 
sions should not be strongly affected if the gas is far from ionization 
equilibrium. It is likely that our spectral extraction region contains 
gas with a range of temperatures and ionization states; however, 
any effect this has on our conclusions below about the gas density 
and shock dynamics is likely to be smaller than the uncertainties in- 
troduced by our geometric assumptions in Section l4~2l In addition, 



the electron and pr oton fluids in the shocked plasma may not be 
in equilibrium (e.g. lKraft et alj |2007a). We discuss this possibility 
further in Section |4~3] 

As we are unable to constrain these effects with our data, we 
adopt the simple, fixed abundance, single-temperature model with 
Gaussian component to account for the prominent Mg line, as dis- 
cussed above. We adopt an abundance of Z — 0.5 times solar to 
compute gas density and pressure, rather than the lower value from 
the free abundance fit, as this is likely to be artificially low due to 
using a single temperature model, as discussed above. As the X- 
ray emissivity is due to both continuum mainly from hydrogen and 
helium and line emission from heavy elements and varies with the 
square of density, the inferred gas density varies at most as Z 1 ^ 2 . 
Thus, the assumed abundance will introduce uncertainties of the 
order of ten per cent at most. 

3.3 Spectral structure of the shell 

We divided the entire lobe into small regions to map the spectral 
structure across the shell. The regions were defined by requiring a 
minimum of 2000 0.4 - 2.5 keV counts (before background sub- 
traction) to allow good constraints to be obtained on the spectral 
parameters. Background spectra were extracted from a region ad- 
jacent to the shell at the same radius as the source region, in order 
to subtract correctly the emission from the galaxy ISM and from 
the wings of the nuclear PSF. This resulted in a total of 50 re- 
gions, which were each fitted separately with APEC (Model III) 
and power-law (Model IV) models in the same way as for the fits 
for larger regions described above. Fig. [5] shows a spectral index 
map for the power-law fits (top panel), and a temperature map ob- 
tained using the APEC fits (middle panel). For each spectral map, 
a corresponding map of the fit statistic (x 2 /n) is also shown. 

As expected from our analysis above, the inner regions of the 
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Table 1. Model fits to X-ray spectra. All fits are joint fits to spectra from all seven observations in the energy band 0.4 - 3.0 keV. Galactic absorption of 
8.4 X 10 20 cm -2 was assumed for all of the fits. For Model V we list the inferred, unabsorbed 0.4 - 3.0 keV fluxes for each model component as a measure 
of the relative importance of thermal and non-thermal emission in each region. 



Model 




Parameter 




Region 1 


Region 2 


Region 3 


Model I: 


APEC (free abundance) 


kT (keV) 






2.03«- 24 4 


0-91±g-gJ 






Z(Z & ) 




< 0.007 


< 0.05 








X 2 (d.o.f.) 




1010 (743) 


344 (308) 


137 (97) 


Model If: 


APEC (fixed abundance) 


kT (keV) 




9 99 +0.08 
z-z -0.06 


9 5q +o.i3 


90+ 03 






Z(Z Q ) 




0.15 


0.15 


0.15 






X 2 (d.o.f.) 




1213 (744) 


377 (309) 


138 (98) 


Model m: 


APEC (fixed abundance) 


kT (keV) 




3 26+ - 09 


3 58+ - 24 


Q7+ 01 

u - a '-0.02 






Z(Z & ) 




0.5 


0.5 


0.5 






X 2 (d.o.f.) 




1478 (744) 


437 (309) 


177 (98) 


Model IV: 


power law 


r 




9 9n +0.02 
Z ' ZU -0.02 


2 01+ 04 


9 94+O.O5 
-0.05 






X 2 (d.o.f.) 




931 (744) 


333 (309) 


445 (98) 


Model V: 


APEC + power law 


kT (keV) 




23+ - 09 


0.93° 


n sq+ 05 

°- 89 -0.06 






Z(Z Q ) 




0.5 


0.5 


0.5 






r 




9 ic+0.06 
z - lo -0.05 


2 00+ 15 
z ' uu -0.06 


9 ,0+0.30 
z - lo -0.27 






X 2 (d.o.f.) 




928 (742) 


332 (307) 


141 (96) 






Fapec 




) x IO" 14 


< 4.8 x 10~ 14 


(2.0 ±0.2) x 10" 14 








(4.7±° i 


) x IO" 13 


(1.38too?) x 10-13 


(9.2+^) x IO" 15 



a The temperature was completely unconstrained for this model. 



shell (close to the AGN) have very different spectral properties to 
those in the outer parts (e.g. Region 1) - for the power-law fits, the 
inner parts have much steeper spectral indices (and often very poor 
fit statistics), consistent with a thermal interpretation for these re- 
gions. In the region between ~ 2 arcmin and ~ 3 arcmin from the 
nucleus none of the fits is good. This is partly because the emission 
there has lower surface brightness than elsewhere, but mainly due 
to contamination from th e wings of the PSF from the bright tran- 
sient source at this radius ( Sivakoff et al .120081) . The outer parts are 
best fitted by a power law, with the spectral index varying from 
~ 1.7 to ~ 2.4, consistent with a predominantly synchrotron ori- 
gin. 

The bottom left-hand panel of Fig. [5] shows a plot of spec- 
tral index as a function of distance from the central AGN. For this 
analysis we defined a set of annuli centered on the nucleus, and ex- 
tracted spectra from regions consisting of the intersection of each 
annulus with a polygon enclosing the shell region. We find that the 
spectral index is steepest in the central regions, where the power- 
law fit is not a good representation of the data. In the outer regions 
the spectral index remains fairly close to a value of 2, but there is a 
statistically significant decline in the spectral index in the last three 
bins. This trend is present whether the results of Model IV or V are 
adopted. This is discussed further in Section l4~2l 

We also investigated the contributions of thermal and non- 
thermal emission to the observed flux as a function of distance 
from the nucleus by fitting Model V (APEC + power law) to the 
same regions as for the spectral index plot discussed in the previous 
paragraph. The bottom right hand panel of Fig. [5] shows the frac- 
tion of the total unabsorbed 0.4 - 3.0 keV flux in the thermal and 
non-thermal model components as a function of distance from the 
nucleus, indicating that the X-ray flux is dominated by non-thermal 
emission except in the very inner regions. The strong evidence for 
dominance by thermal emission in Region 2 leads us to conclude 
that this is likely to be a real difference in the shell emission, rather 
than an artefact of poor background subtraction. 



4 DISCUSSION 

Our spectral analysis of the X-ray shell emission strongly favours 
a non-thermal interpretation for the majority of the X-ray emis- 
sion in the outer half of the shell, due to the unrealistically low 
abundances required by thermal models, as well as the lack of ex- 
pected emission-line signatures, particularly the Fe L shell com- 
plex. These results, together with the strong resemblance between 
the highly edge-brightened shell emission and the X-ray emission 
fro m synchrotron-dominate d supernova remnants, such as SN1006 
(e.g. lRothenflug e t al. 2004), lead us to conclude that an X-ray syn- 
chrotron origin for the majority of the X-ray shell emission from 
Cen A is viable. Further constraints on this model come from lim- 
its at other wavelengths. It might be expected that a plausible en- 
ergy distribution for an electron population producing X-ray syn- 
chrotron emission would also produce detectable emission at ra- 
dio, optical, UV or infrared wavelengths, and so the model must 
pass this test. It must also be possible to construct a self-consistent 
picture of the radio-lobe dynamics, which explains both the strong 
thermal emission from Region 3 and the lack of any thermal sig- 
natures in the outer rim region. Below we discuss the constraints 
on a synchrotron model that can be obtained from multiwavelength 
measurements and from an analysis of the shock conditions and 
lobe dynamics. We also consider the dependence of our conclu- 
sions on the assumption of electron-ion equilibrium. Finally, we 
consider the implications of our results for particle acceleration in 
Cen A and other radio galaxies, and for Cen A as a source of ultra- 
high energy cosmic rays (UHECRs) and TeV emission. 



4.1 The broad-band SED of the shell region 

The brightest part of the X-ray shell surrounding Cen A's SW lobe 
is the outer rim, which is spatially separated from the region coinci- 
dent with the radio emission from the radio lobe itself. If the X-ray 
emission from the shell is produced via the synchrotron process, 
then it would be expected to originate in an electron population 
that extends to energies much lower than those emitting in the X- 



Cen A radio-lobe shock 7 



St 






Distance (arcsec) Distance (arcsec) 

Figure 5. Spectral structure across the shell, as discussed in Section 13.31 Top: spectral index map, assuming Model IV (power law) for all regions (left: 
blue=1.5, red=4.0), with corresponding map of x 2 / n (right: red = 0.5, blue=2.5). Middle: temperature map, assuming Model III (APEC with Z = 0.5) for all 
regions (1: red=4.0 keV, yellow=0.7 keV), with corresponding map of x 2 / n ( r: re ^ = " 5, blue=2.5). Bottom left: spectral index as a function of distance from 
the nucleus. The large error bar for the third data point is because a large fraction of this region is excluded due to a bright X-ray transient present in several 
datasets (Sivakoff et al. 2008). Bottom right: fraction of flux in thermal (red) and non-thermal (blue) components for Model V as a function of distance from 
the nucleus. 



ray band. No shell emission has been reported previously at radio, 
optical, IR, UV or radio wavelengths. We determined limits in the 
radio, infrared and UV using VLA, Spitzer and GALEX observa- 
tions to construct a broadband SED (Fig. |6j for a region of shell 
emission surrounding the radio lobe (we used Region 2, discussed 
above, which is exterior to the radio lobe in projection, so as to 
obtain as tight a radio constraint as possible). 

The radio limit shown in Fig. [6] was obtained from the 1.5- 



GHz radio map described in Section [2] and is a 3cr limit based on 
the off-source rms noise, as there is no detection of radio emission 
from Region 2. We used the S pitzer IRAC and MIPS dBrookes et al. 
l2006l :[Hardcastle et al. 2O06) and GALEX observations jNeff et al. 
120031 : Hardcastle et al. 2006) to obtain infrared limits at 3.6, 4.5, 
5.8, 8.0 and 24 /im and a UV limit at 231 nm. Th e analysis was 
carried out as described in lHardcastle etal] d2006h . Limits were 
calculated for Region 2, using a larger region immediately outside 
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Figure 6. The broadband SED for Region 2, including upper limits at 1.4-GHz, 3.6, 4.5, 5.8, 8.0 and 24 fim, 231 nm, and the measured X-ray flux density at 1 
keV and spectral index. Solid line is a synchrotron model fit with (1) a = 0.5, 7 m i n = 10, "/break = 4.4 X 10 4 , 7 moI = 5 X 10 9 (discussed in Section BTt . 
and (right panel) a = 0.5, no break, and "/max = 3 X 10 s (discussed in Section l4^4V 



Region 2 as a background. At all bands from the IR to the UV the 
dominant background is stars, and so this choice of background 
will give a strict upper limit. As an alternative strategy, we also 
used galaxy-subtracted IR images to obtain tighter limits on the 
IR emission from the shell. This allowed us to decrease the upper 
limit at 3.6 /im by roughly an order of magnitude. Similar improve- 
ments could be obtained using this method at other wavelengths, 
but would not significantly affect the constraints on a synchrotron 
model, as discussed below. Figure [6] includes all of the limits ob- 
tained, including the galaxy-subtracted limit at 3.6 /im. 

Fig. [6] shows a synchrotron model chosen to fit through both 
the X-ray data point and spectral index and the 1.4-GHz upper 
limit. The plotted model assumes a = 0.5, jmin = 10, "/break = 
4.4 x 10 4 , -fmax = 5 x 10 9 , and uses the 1.4-GHz upper limit 
to normalise the electron spectrum, assuming equipartition in the 
shell material. The volume enclosed in Region 2 was modelled by 
approximating to a box of dimensions 573 x 7 x 46 arcsec, where 
the length and width are the dimensions of the extraction region, 
and the depth is approximated by that halfway between the inner 
and outer curved boundaries of the regions, assuming a spherical 
shell of appropriate radius to fit the shell curvature. This gives a 
volume of ~ 3 x 10 64 cm 3 . None of the infrared or UV limits can 
rule out the synchrotron model described above. In addition, since 
the radio data point is an upper limit, a model with a significantly 
higher break frequency (as typically fo und for supernova remnants 
with X-ray synchrotron emission, e.g. iDver et al]|2004l ; iRho et all 
120031 ; see also Section |4~4t , which could be orders of magnitude 
below the IR limits, also cannot be ruled out (see the right hand 
panel of Fig. [6j. Hence we conclude that the broadband SED for 
the rim of the Cen A SW shell is consistent with a synchrotron 
interpretation for the X-ray emission. 

Although a non-thermal contribution to the spectrum in the 
NE part of the shell (Region 3) is not required, the limits on 
a power-law contribution to the Region 3 spectrum from fitting 
Model V are consistent with a significant flux from synchrotron 
emission in this region (the fit for Model V implies an unabsorbed 
1-keV flux density from the power-law component of ~ 2 njy). 



Therefore, we cannot rule out significant particle acceleration in 
Region 3 as well. 

Finally we also investigated whether inverse Compton emis- 
sion could explain some or all of the shell X-ray emission. X-ray 
emission from IC scatting could be significant in the regions coin- 
cident with the radio lobes (whose X-ray emission we have hitherto 
been interpreting as coming from the shell material in front of and 
behind the lobe). We therefore modelled the electron distribution 
in the radio lobe in order to predict the expected level of IC emis- 
sion. We approximated the lobe as a cylinder of length 270 arcsec 
and radius 80 arcsec, using radio flux measurements at 240 MHz 
and 8 GHz to determine the shape and normalisation of the elec- 
tron distribution, assuming equipartition. We assumed "/ m in ~ 10, 
as above, and chose "/max to fit the observed radio spectrum, giv- 
ing a value of 1.2 x 10 4 . The predicted 1-keV flux density for the 
lobe from inverse-Compton scattering of the cosmic microwave 
background radiation (IC/CMB) (the dominant process) and syn- 
chrotron self-Compton (SSC) is ~ 3.3 njy, which is factor ~ 30 
below the observed flux from the lobe region. Hence it is possible 
that a small fraction of the observed X-ray flux from Region 1 is 
produced by the IC/CMB process. For the lobe emission to be pro- 
duced entirely by IC/CMB would require a magnetic field a factor 
~ 6 below equipartition. An inverse-Compton model for the edge- 
brightened rim of the X-ray structure (e.g. Region 2) is implausible, 
due to our tight limits on the radio emission in this region, and be- 
cause none of the m any examples of IC/CMB-detected lobes (e.g. 
ICroston et ai]|2005l) show edge-brightening. We therefore consider 
it unlikely that our conclusions are affected by possible contamina- 
tion from X-ray IC emission. 

4.2 Lobe expansion and shock dynamics in a synchrotron 
model 

As the X-ray emission from the north-east region of the shell (Re- 
gion 3) is unambiguously thermal in origin (see Section [3~2t , this is 
the best region to use for investigating the shock jump conditions. 
We approximated the volume of shell enclosed within Region 3 as a 
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Table 2. Adopted physical parameters for Regions 1 and 3 



Parameter 


Region 2 


Region 3 


T 3hM (keV) 


11.0 


0.95 


T ism (keV) 


0.35 


0.35 


Z (Z ) 


0.5 


0.5 


n p ,shell (cm -3 ) 


0.004 


0.033 


Hp, ism (cm ) 


0.001 


0.01 


Pshell (dyne cm -2 ) 


1.1 x 10" 10 


1.1 x lO" 10 


Pism (dyne cm -2 ) 


1.3 x 10~ 12 


1.1 x 10- 11 



box of dimensions 0.480 x 0.127 x 0.769 arcmin, where the length 
and width are the dimensions of the spectral extraction region, and 
the depth along the line of sight was taken to be the chord of a circle 
having the curvature of the outer edge of the shell in the north-east 
region at a distance halfway between the inner and outer edges of 
the spectral extraction region. This leads to an assumed volume 
of 1.7 x 10 63 cm 3 . With this volume, and assuming the APEC + 
Gaussian model described in Section [3T2l we obtain a density of 
n p — 0.033 cm -3 , which corresponds to a pressure in the shocked 
gas of 1.1 x 10 -10 dyne cm" 2 , assuming kT — 0.95 keV as dis- 
cussed above. The shell density in Region 3 is roughly 3.5 times 
the density of the ISM at this distance, based on the beta mode l 
fitted to the external gas density distribution bv lKraft et al.l d2003h . 
although we note that the ISM density is not well constrained, be- 
cause the distribution of gas is likely to deviate significantly from 
spherical symmetry, and to be inhomogeneous in the inner regions 
of the galaxy. Hence the density constrast is consistent with the 
value of 4 expected for a s trong shock, as was concluded previ- 
ously bv lKraft et all d2007ah . If the hot-gas shell and radio lobe are 
assumed to be in pressure balance, then the pressure inside the ra- 
dio lobe is expected to be ~ 1.1 x 10~ 10 dyne cm~ 2 , as calculated 
above, which is about an order of magnitude higher than the lobe 
minimum pressure. This is of the same order as the differences be- 
tween lobe minimum pre ssures and external pressu res for FRI radio 
galaxies in general (e.g. ICroston et al.l |2007| . [2008 ) . The pressure 
contrast between the shell and the ISM is P s hm /Pism = 10, us- 
ing th e gas density distribution and temperature from iKraft et al.l 
(2003). Assuming ram-pressure balance and that 7 = 5/3, this 
pressure jump implies that, at its north-east edge, the radio lobe is 
expanding at a Mach number of ~ 2.8, which corresponds to an 
exp ansion speed of ~ 850 km s . 

IKraft et all J2003I) found that the apparent density contrast be- 
tween shocked and unshocked gas in the outer parts of the shell 
(assuming a thermal interpretation for all the shell X-ray emission) 
is much higher than the prediction for a strong shock, and derived 
the expansion speed assuming that the observed shell is shocked 
gas that has cooled rather than the component that directly enters 
the jump equations for shocked gas. A non-thermal model for most 
of the X-ray emission within the shell can provide a simpler in- 
terpretation of the larger X-ray-intensity constrast inside and out- 
side the shell. We checked the low apparent contribution from ther- 
mal emission in the outer regions of the shell is consistent with 
the thermal emission expected from a shell of shocked gas a fac- 
tor of 4 denser than the ISM in the shell's outer regions (for sev- 
eral assumed temperatures). We assumed an unshocked ISM den- 
sity of 1.0 x 10~ 3 cm -3 at the outer edge of the shell (Kraft et 
al. in prep.: note tha t thi s is a factor 1.7 lo wer than that assumed 
bv lKraftetal]|2003l and lKraft et al]|2007bl) . so that the shell gas 
would be expected to have a density of ~ 4 x 10 -3 cm -3 . We 
calculated the volume corresponding to the outer shell region dis- 



cussed in Section [3~T1 assuming a half spherical shell of outer ra- 
dius 121 arcsec. We estimated the thickness of the shell using the 
1.4-GHz radio map and the merged 0.4-2.5 keV event file by mea- 
suring the separation between the edge of the radio lobe and the 
brightest part of the shell edge. We took the average of three mea- 
surements around the edge of the shell, giving a mean thickness of 
16.3 ± 1.5 arcsec, where the measurement uncertainty (which may 
be partly due to real variations in the thickness of the shock front) 
is larger than any astrometric uncertainty or resolution effects. Our 
measured s hell thickness is c onsistent with the value of 17 arcsec 
adopted bv lKraft et all J2007a). This assumed geometry leads to an 
expected emission measure of ~ 1.7 x 10 58 cm -3 (corresponding 
to a 0.4 - 3.0 keV flux of 1.3 x 10~ 14 erg cm -2 s" 1 ). The measured 
normalization for the APEC component of Model V discussed in 
Section [3~T| corresponds to an emission measure more than thirty 
times higher, hence the data are consistent with this emission being 
present. However, as the APEC component in the spectrum of this 
region is not required by the data, the spectra can also be well fit- 
ted by an APEC + power-law model with the APEC normalization 
fixed at the value corresponding to the expected emission measure, 
for any reasonable temperature (e.g. between kT = 0.5 and 15.0 
keV - our adopted shock model as discussed below implies a post- 
shock temperature of ~ 11 keV). In all cases the thermal contri- 
bution to the fit is negligible. Hence the X-ray spectra for the outer 
shell are consistent with a model in which a strong shock is occur- 
ring, with the X-ray emission dominated by synchrotron emission 
from high-energy particle acceleration. 

If we assume that the radio lobe is isobaric, and that the ra- 
dio lobe and shocked gas are in pressure balance everywhere, then 
the pressure in the shocked gas shell at its outer edge (e.g. Region 
2) should be the same as the lobe internal pressure and shell pres- 
sure in Region 3, which we calculated above. This implies that the 
expansion speed at the outer edge of the lobe where we see the 
brightest non-thermal emission should be much higher, given the 
significantly lower external gas density (~ 10~ 3 cm' 3 ). Assuming 
an ISM gas temperature of 0.35 keV (Kraft et al. in prep), we ob- 
tain a pressure jump P s hM /Pism ~ 87, corresponding to a Mach 
number of ~ 8.4, and a lobe expansion speed at its o uter edge of 
~ 260 km s _1 . This speed is similar to that inferred bv lKraft et al.l 
( 2003). We can also infer a transverse expansion speed for the radio 
lobe similarly, us ing the external d ensity at the lobe midpoint given 
by the profile in lKraft et alj (2003D . which we find to be ~ 1400 km 
s _1 . The observed axial ratio of the radio lobe is roughly consis- 
tent with the inferred longitudinal and transverse expansion speeds 
for the lobe if the expansion has been self similar throughout the 
source's lifetime (though in practice this is unlikely to have been 
the case). The adopted physical parameters for Regions 2 and 3 are 
summarized in Table [2] 

As for Region 3, we have assumed here that 7 = 5/3. This 
will be correct if the shell pressure is dominated by the thermal 
component of the shocked gas, rather than by the pressure in the 
relativistic particles that are producing the observed X-ray syn- 
chrotron emission. As we concluded above, the observations are 
consistent with a scenario in which the hot shocked gas has the den- 
sity expected for a strong shock, but is not separable from the domi- 
nant emission attributed to synchrotron radiation. However, we also 
checked whether the internal pressure of the synchrotron-emitting 
particles is likely to significantly alter the shock dynamics. Assum- 
ing the SED discussed in Section l4~Tl with an equipartition magnetic 
field of ~ 8 /iG and no energetically significant protons, we find 
that the pressure in cosmic rays and magnetic field is ~ 2 x 10 -12 
dyne cm -2 , which is insignificant compared to the thermal pres- 
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sure. We note that this is an upper limit on the pressure in the rela- 
tivistic particles if in equipartition, because the SED model is based 
on upper limits at radio, IR and UV wavelengths. Magnetic field 
strengths of ~ or higher, more than an order of magnitude 

higher than the equipartition value, would be required for the non- 
thermal (magnetic) pressure to dominate, assuming j m in = 10. A 
departure from equipartition in the direction of electron dominance 
could also lead to non-thermal pressure dominating over thermal 
pressure; however, this would also have to be extreme. 

We can also estimate the total energy of the radio source and 
the power of the jet. Approximating the shock front as a prolate 
spheroid with semi-major axis 167.2 arcsec and semi-minor axis 
115.0 arcsec, we find an enclosed volume of 1.6 x 10 66 cm 3 . The 
cavity enthalpy, APV is then 7 x 10 56 ergs. We make a rough es- 
timate of the kinetic energy of the shell using the expansion speed 
midway along the shell, and the gas density for Region 1, to ob- 
tain Ek = 1.2 x 10 56 ergs. Assuming a constant expansion speed, 
we obtain a source age estimate of ~ 2 x 10 6 y, and hence we 
estimate a mean jet power over the source lifetime of ~ 10 43 
erg s _1 . We can also estimate an instantaneous jet power as PV, 
where V = V(a/a + 26/6), with a and 6 the semi-major and 
semi-minor axes respectively. If the lobe expansion is self-similar, 
then V — ZVa/a — Anb 2 a. As one end of the expanding lobe is 
tied to the AGN, 2d = v s hock, the shock expansion speed at its 
outer edge, and hence the instantaneous jet power, p, is given by 
p — 2TvPb 2 v s hock = 6.6 x 10 42 erg s _1 (using the shock speed 
of 2,600 km s _1 at the outer edge of the shock). This is roughly a 
factor 15 times hi gher than the current radiative power of the AGN 
jEvans et alj2004h . 

4.3 Electron-ion equilibrium in the shocked gas 

Our conclusions about the shock dynamics in Section |4~2l depend 
strongly on our adopted spectral model for the thermal emission in 
Region 3. As discussed in Section [3~2l the region shows a some- 
what unusual spectrum, with a prominent magnesium line, that is 
not well fitted by a single APEC model. As the gas in this region 
is assumed to have passed through the shock front recently, it is 
likely that the ionization balance has been affected by the shock, as 
discussed in Section [3~2l It is also likely that our spectral extrac- 
tion region contains gas with a range of temperatures and ioniza- 
tion states, but in Section [3~2l we concluded that any effect on our 
conclusions about the gas density and shock dynamics is likely to 
be smaller than the uncertainties introduced by our geometric as- 
sumptions in Section |4~2l However, an additional complication is 
that the electron and ion fluids may not have achieved equilibrium, 
as was discussed in th e context of earlier work on the Cen A shock 
bv lKraft et all d2007 al) . 

For gas temperatures of ~ 1 keV and below, and photon ener- 
gies below ~ 1 keV, at the energy resolution of ACIS-I the thermal 
continuum is submerged by lines, particularly those of the Fe L 
complex. This makes fitted temperatures insensitive to the contin- 
uum. If the gas is not in ionization equilibrium, the relationship 
between line strengths and gas temperature is model-dependent, so 
that a detailed model of the shock and subsequent return to ion- 
ization equilibrium is required to determine the gas temperature 
accurately. On the other hand, the shape of the continuum is deter- 
mined largely by the electron temperature. In the absence of a de- 
tailed model for the ionization state of the gas, we have attempted 
to constrain the postshock electron temperature by fitting the spec- 
trum of Region 3 in the energy range 1.5-3 keV, where the lines 
are weaker relative to the continuum. A 1.3 keV Gaussian line was 



included in this fit, as discussed in Section [3"!2l Its energy and nor- 
malization were fixed by fitting it together with an APEC model 
over the energy range 0.5 - 2 keV. For the 1.5-3 keV fit, the abun- 
dances were fixed at 0.5 solar. The resulting 90% confidence range 
for the temperature is 0.52 - 1.39 keV, consistent with the previous 
thermal fit. 

If the flow time through Region 3 is comparable to or smaller 
than the electron-ion equilibration time, then the shoc ked protons 
there c ould be significantly hotter than the electrons (cf. lKraft et al.l 
l2007al) . In that case, the postshock pressure may be significantly 
greater than estimated above, only providing a lower limit on the 
true pressure of the radio lobe. The value of the lobe pressure is 
important, since it determines the strength and speed of the shock, 
and the estimated jet power scales approximately as P 3 ^ 2 , so this 
issue warrants careful scrutiny. 

To determine how much the electron and proton tempera- 
tu res may differ, we h ave used the two-fluid hydrodynamic model 
of iKraft et al.l ( 12007 al) . Motivated by observations of Galactic su- 
pernova remnants, this model assumes that the protons are heated 
rapidly in a collisionless shock, while the electrons are heated by 
Coulomb collisions with the protons. The model is spherically sym- 
metric and, for the cases illustrated here, the shock is driven by 
constant power input into an expanding "radio lobe" at the centre 
of an initially uniform atmosphere. For scaling purposes, the shock 
radius is taken to be 1.6 kpc (90 arcsec). 

In terms of the electron and proton temperatures, T and T p , 
respectively, the energy exchange can be expressed as dT a /dt — 
v(T p — T c ), where the rate coefficient is determined chiefly by 

3/2 

the electron temp erature, approximately v oc n c T a (e.g. 
IKraft et ai]|2007al) . Although the electron-ion equilibration time is 
relatively short, ~ 2 x 10 5 y for n p = 0.033 cm" 3 and kT = 0.97 
keV, it is not much shorter than the sound crossing time of Region 
3, ~ 2.6 x 10 J y. Increasing the temperature of the shocked gas 
increases the electron-proton equilibration time while also decreas- 
ing the sound crossing time. Thus, although electrons and protons 
in Region 3 should have much the same temperature if the electron 
temperature there is close to 1 keV, their temperatures can differ 
significantly if the electron temperature is moderately higher. 

This is illustrated in Fig. [7] which shows postshock proton 
and electron temperature profiles for two cases. In the first case 
(Fig. |7Ji), the emission measure weighted electron temperature in 
Region 3 (1.48 - 1.60 kpc) is 1 keV and the RMS proton density is 
0.033 cm~ 3 . In this case the difference between the electron and 
proton temperatures is negligible. For the second case (Fig. |7p), 
the RMS density is the same, but the emission measured weighted 
electron temperature is close to 1.4 keV, the 90 per cent level up- 
per limit. Here, the proton temperature is ~ 70 per cent greater 
than the electron temperature. The postshock pressure would then 
be roughly 35 per cent larger than suggested by the electron tem- 
perature. This disparity between electron and proton temperatures 
grows rapidly with electron temperature. However, if the electron 
temperature is smaller than the 1.4 keV upper limit, our estimate 
for the postshock pressure should be a reasonable estimate of the 
true pressure. 



4.4 Implications of the X-ray synchrotron interpretation for 
particle acceleration at radio-lobe shocks 

In Section |4~T| we demonstrated that a synchrotron origin for the 
shell X-ray emission is consistent with the broad-band SED of the 
shell emission in a region exterior to the radio lobe (where the most 
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Figure 7. Post-shock proton temperature (solid line) and electron temperature (dotted line) profiles obtained from our two-fluid hydrodynamical modelling 
(Section|43) for kT = 1 keV (left) and fcT = 1.4 keV (right). 



constraining radio limits could be obtained). We have also shown 
that the lack of thermal X-ray emission from the outer parts of the 
shell is consistent with the expected level of emission from a strong 
shock - the observed thermal X-ray emission in the more cen- 
tral NW shell region (Section [3.2t and the lack of any signatures 
of such emission in its outer parts (Section |3. It are as expected 
for a strong shock propagating into Cen A's environment. We con- 
clude that the majority of the shell X-ray emission is synchrotron 
emission from high-energy particle acceleration at the shock front. 
Interestingly, the inferred expansion speed of the radio lobe at its 
outer edge, as calculated in the previous section, is very similar 
to the expansion speeds of supern ova remnants (SNRs) that show 



X-ray synchrotron emission ( e.g. IVink et al.l 120061 ; I Warre n et al.\ 
12003 ; lRothenflugetai]|2004 lAcero et alj|2007l) . whereas at the 
NW edge, where no significant non-thermal emission is detected 
(although particle acceleration to X-ray-emitting energies cannot 
be ruled out - see Section[4j}, the expansion speed is significantly 
lower. If the inner lobes of Cen A are indeed expanding at these 
high velocities, then it is perhaps not unexpected to see such a 
strong signature of high-energy particle acceleration. 

The statistically significant decrease in spectral index in the 
outer ~ 50 arcsec of the shell (Section [3.3b could indicate differ- 
ences in the particle acceleration behaviour as a function of distance 
and hence of shock expansion velocity. A possible interpretation 
is that the spectrum is flattest at the outermost edge of the shell 
(where the expansion speed is highest) because of a higher value 
for 7 ma a;, the observed high-energy cut-off of the electron popula- 
tion, with the high-energy cut-off moving to lower energies, thus 
steepening the spectrum, in regions closer to the nucleus where the 
shock front is slower due to a denser environment. This would re- 
quire 7 mal in the range 10 7 — 10 9 over the bright shell region, and 
significantly lower than 7 x 10 7 in Region 3. (We note that, assum- 
ing the observed value of "/max is determined by radiative losses, 
the high-energy cut-off for hadrons would be significantly higher). 
If the X-ray-emitting electrons are being accelerated by diffusive 
shock acceleration, and if the high-energy cut-off is determined by 
the synchrotron loss timescale, then we can determine the depen- 
dence of 'ymax on shock speed, as well as the expected value of 
7 mol for the shock conditions we observe. The high-energy cut-off 
for particle acceleration has been discussed extensively in the con- 
text of X-ray synchrotron emission from supernova remnants (e.g. 
lLagage & Cesarskvll983l ; lRevnoldsll 996). For electrons, the accel- 
eration timescale to a given energy E depends on Bi (the pre-shock 



magnetic field strength), E, v a hock, the ratio of the electron mean 
free path along the magnetic field to the gyroradius (/, assumed to 
between 1 and ~ 10), and on a factor Rj that takes into account the 
decrease in efficiency of scattering if the magneti c field direction i s 
not parallel to the shock. We used equation (1) of lRevnoldsl dl996h : 



E„ 



0.1(/i?jBi)- 1/2 w 8 ergs 



(1) 



where units are cgs and the shock speed is in units of 10 8 cm s" 1 . 
We assumed / = 10, Rj = 1 and Bi = 1/iG, and find that 
Jmax ~ 10 8 for the outer edge of the shell and ~ 3 x 10 7 at Region 
3. In order to fit a synchrotron model consistent with the observed 
X-ray spectral index for Region 2, and with equipartition B fields, 
we require ~fmax in the range 2 x 10 8 — 4 x 10 8 (see Fig.|5J, which 
is roughly consistent with this analysis. This suggests that a model 
in which the observed X-ray spectral changes can be explained by 
changes in ymax as a function of shock speed is viable. We have 
not considered the role of changes in B\ and / along the shell in 
this discussion - these are likely to vary somewhat with position as 
well, but the dependence of ^max on these parameters is weaker so 
that the change is shock speed is likely to be the dominant effect. 



4.5 Implications for large-scale shocks in other galaxies 

If we accept a synchrotron origin for the X-ray emission from Cen 
A, then the X-ray properties of the outer shell region cannot be 
used to measure the impact of the radio galaxy on its environ- 
ment directly; however, the detection of thermal emission from the 
shocked gas in Region 2 allows us to infer the dynamics of the 
shocked shell and radio lobe indirectly as discussed above. The 
results on the overall dynamics of the radio lobe and its environ- 
m ental impac t are n ot significantly different from the conclusions 
of lKraftetai] ( l2003h . Unlike the shock in Cen A, the strong sho ck 
detected in the similar system NGC3801 dCroston et"al]|2007l) is 
clearly dominated by thermal X-ray emission, with no significant 
X-ray synchrotron emission. It is likely that the lobes of NGC 3801 
are expanding much more slowly than the outer parts of Cen A 
(~ 750 km s _1 , similar to Cen A's expansion speed at Region 3), 
and so the lack of strong high-energy particle acceleration in that 
system is probably explained by this lower speed. 

The bolometric X-ray luminosity of the shell emission in Cen 
A is ~ 4 x 10 39 erg s _1 , and so shocks on similar physical scales 
will be undetectable in more distant radio galaxies; however, these 
results could also have implications for powerful FRII radio galax- 
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ies, whose lobes are much larger. A number of nearby FRII radio 
galaxies have be en found to resid e in com paratively gas poor en- 
vironments (e.g. lKraftetalJl200ll2007bh . and so may be having 
similarly dramatic effects on their environments. If supersonic ex- 
pansion is occurring in FRIIs at speeds similar to that of Cen A, 
then high-energy particle acceleration could be important in FRII 
environments. We considered the example of the nearb y FRII radio 
galaxy 3C 33, whose poor environment is discussed in lKraft et al.l 
( I2007bf) . If the lobes of 3C 33 are expanding at a similar speed to 
those of Cen A, and we assume (i) that the particle acceleration 
is similarly efficient and (ii) that the density of accelerated parti- 
cles is proportional to the density of the shocked gas, then the ex- 
pected X-ray synchrotron flux should scale as S oc riismVDj 2 . 
The external density ratio 713033 /ncenA is ~ 0.01, the volume ra- 
tio is ~ 7500, and the square of the ratio of luminosity distances is 
1.9 x 10 -4 , giving an expected reduction in X-ray flux of a factor 
~ 70 and a predicted 1-keV flux density from synchrotron emis- 
sion of ~ 1 njy. This is comparable to the level of non-thermal 
X-ray emission measured by Kraft et al. and attributed to inverse- 
Comp ton scattering of the cosmic microwave background (CMB) 
(e.g. lCroston et alj2005h , and so for 3C 33 the observed X-ray flux 
levels from the radio lobes do not rule out the possibility that high- 
energy particle acceleration is occurring at a similar level to Cen 
A. However, the spectral index for the non-thermal X-ray emission 
from 3C 33 and other FRII radio lobes (~ 1.5 — 1.7) is flatter than 
that of the shell in Cen A, and no edge-brigh tening is observed for 
the X-ray lobes of FRIIs dCroston et al]|2005h . In addition, the X- 
ray flux from this process decreases as (1 + z) 4 , unlike the IC/CMB 
flux, and so IC/CMB is predicted to dominate by a large factor in 
the majority of FRIIs, which are more distant than 3C 33. Never- 
theless, it is plausible that particle acceleration by expanding FRII 
lobes may result in significant cosmic ray populations and poten- 
tially detectable X-ray synchrotron emission for sources in weaker 
environments. 



4.6 Implications for Cen A as a cosmic ray source 

It is interesting to consider the implications of these results for our 
understanding of Cen A as a source of cosmic rays, particularly in 
light of recent results from the Pierre Auger Observatory showing 
an overdensity of ultra-hig h energy cosmic rays (UHECRs) from 
the d irection of Cen A dAbraham et al\ 120071 : [Moskalenk o et al.l 
2008). To have efficient cosmic-ray acceleration to high energies, 
the particle gyroradius cannot be significantl y large r than the size 
of the acceleration region. By the criterion of Hillas (1984), for an 
assumed acceleration region of ~ 300 pc (the observed thickness 
of the shock front: Section l4~2l > and B = 7/iG, typical particle en- 
ergies cannot exceed ~2x 10 18 eV. Hence magnetic fields signifi- 
cantly higher than equipartition (e.g. ~ 400/iG) would be required 
for the shock front to accelerate cosmic rays to 10 20 eV. There 
is considerable evidence that magnetic fields ca n be amplified by 
large factors at the f orward shocks of SNRs (e.g. lBerezhko & Volkl 
2004; Ellision et alj|2004l) (for example lower limits on B in Cass 
A a re of the order of thos e required here to confine CRs of 10 2 ° 
eV: IVink& Laming! 120031) . and so if similar processes were op- 
erating at the Cen A shock front then UHECR energies could be 
achieved. However, the equipartition magnetic field strengths in 
SNRs are comparable to values inferred from the thickness of the 
shock front. This suggests that magnetic field amplification by sim- 
ilar factors at the Cen A shock may be unlikely, as it would require 
a large departure from equipartition in the direction of magnetic 
field dominance. The Cen A shock front appears resolved in our 



Chandra observations, with a width of ~ 2 arcsec (correspond- 
ing to ~ 35 pc). We can use the shock thickness to estimate the 
magnetic field strength in the shell, assuming it i s determined by 
synchr otron losses. Following the method used bv lVink & Lamina 
(2003) for the Cas A SNR, assuming the accelerated electrons 
travel away from the shock front with the shocked gas at a speed 
of 0.25v s hock, we find synchrotron lifetimes of ~ 5 x 10 4 y are 
required. These lifetimes for 1-keV photons imply magnetic field 
strengths of ~ 1/iG, consistent with the equipartition values di- 
cussed above, but much lower than would be required to accelerate 
UHECRs. In addition, as discussed above, our X-ray observations 
suggest that we may be observing the high-energy cut-off in the 
electron distribution directly, at energies of ~ 7 = 10 7 - 10 8 , 
which corresponds to electron energies of ~ 10 14 eV. Neverthe- 
less, if the cut-off energy is determined by losses then its value will 
be higher for hadrons. We conclude that it is possible that the Cen A 
shock could accelerate cosmic rays to UHECR energies; however, 
this would require magnetic field amplification by large factors, and 
magnetic domination of the energetics of the synchrotron-emitting 
plasma. The ~ 600-kpc-scale giant outer lobes of Cen A may be a 
more likely candidate for the origin of the UHECRs apparently as- 
socia ted with Cen A (e.g. [Moskalenko et al. 2008; Hardcastl e et al.l 
l2008h . 

4.7 Implications for Cen A as a TeV source 

It is also important to consider the implications of our results for 
emission at higher energies. The detection of X-ray synchrotron 
emission from the shell implies the existence of electrons with 
TeV energies (for plausible values of B), and hence it is possi- 
ble that the shock could produce significant emission at very high 
energies. We therefore investigated the predicted inverse-Compton 
emission from the Cen A shell. The dominant photon fields in the 
shell region that can be scattered to gamma-ray energies are (1) 
starlight from the host galaxy, and (2) the cosmic microwave back- 
ground (CMB). We determined the photon density from starlight 
at the position of th e shell by dep r ojecti ng the V-band surface 
brightness profile of van den Bergrj dl976l) . using the method of 
iMellier & Mathezl d 19871) " to obtain a V-band spectral energy den- 
sity at a radius of 268 arcsec (midway between the inner and 
outer radii encompassing Region 1) of 1.6 x 10~ 28 J Hz" 1 m~ 3 . 
A spectral energy distribution fo r the en tire galaxy was obtained 
using the UV measurements of JWelchl dl979h . optical measure- 
ments from Ide Vaucouleurs dl99ll) . infrared measurements from 
2MASS dJarrett et alj 120031) and IRAS dGolombek etafl 1 19881) , 
which was normalised to the V-band energy density given above. 
We note that the corresponding total energy density at the ra- 
dius of interest is i n good agreement with that determined by 
IStawarz et al. I d2006l) . The analysis was carried out for Region 1, 
so as to include most of the X-ray-emitting electrons, and the elec- 
tron population was modelled as described in Fig. 14.11 We u sed the 
SYNCH code of lHardcastle. Birkinshaw and Worr all ( 1998) to de- 
termine the predicted inverse-Compton emission from starlight and 
IC/CMB. Fig. [8] shows the inverse-Compton predictions for mag- 
netic field strengths of 1 /jG (the value obtained above from consid- 
eration of the width of the X-ray shell) and 7/iG (the equipartition 
value for Region 2 discussed in Section |4~Tt . Also plotted is an up- 
per limit from the High-Energy Stereoscop ic System (HESS) of < 
5.68xl0~ 12 phcirT 2 s~ 4 above 0.19 TeV dAharonian et al.l2005l) . 
which corresponds to a flux density at 0. 19 GeV of 3.8 x 10 -15 Jy 
for an assumed photon index of 2.0, consistent with our assumed 
spectrum. For B = 7/iG, we predict a total flux density at 0.19 
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Figure 8. Predictions for inverse-Compton emission from the SW shell of Cen A from scattering of starlight (orange dash-dotted line) and the CMB (blue 
dotted line). The electron distribution is as in the left-ha nd panel of Fig.jjjQfor B = l^tG (left) and (right). The star indicates the measured X-ray flux 
density, and the upper limit is the HESS measurement of Ahar onian et all (2005). 



GeV from IC scattering of starlight and the CMB of 3.2 x 10~ 16 
Jy, an order of magnitude below the HESS limit, which corresponds 
to F(> 0.19GeV) ~2x 10" 13 photons cm" 2 s" 1 above 0.19 
TeV (for V = 2). However, with a lower magnetic field strength, as 
is suggested by the calculation in Section l4~6l we find a predicted 
total flux density at 0.19 GeV of 1.6 x 10~ 14 Jy, corresponding 
to a total F(> 0.19GeV) ~ 2 x 10" 11 ph cm" 2 s~\ i.e. we 
predict a flux that is comparable to or higher than the HESS up- 
per limit. Hence it is reasonable to think that a TeV detection of 
the Cen A shell may be possible with deeper observations. It is un- 
clear whether the Cen A shell would be expected to dominate over 
emission from the jet (which a lso contain electrons of TeV ener- 
gies: e.g.lHardcastle e t al.l2007h at TeV energies. A full calculation 
of the predicted TeV emission from the jet is beyond the scope of 
this paper; howev er, a rough estimate base d on the X-ray jet prop- 
erties discussed in Hardcast le et all d2007l) suggests a level of TeV 
emission comparable to that from the shell. In principle it would 
be possible to distinguish spectrally between IC emission from the 
jet and the shell, as the diffuse X-ray synchrotron emission in the 
jet has a significantly steeper spectral index than the shell. Hence 
deeper TeV observations of Cen A have the potential to provide 
important constraints on the properties of the fluid in the shocked 
shell as well as the jet. 



5 CONCLUSIONS 

We have presented a detailed study of the X-ray emission associ- 
ated with the south-west inner lobe of Cen A. Our deep observa- 
tions as part of a Chandra Very Large Programme have enabled 
us to establish firmly that the majority of the X-ray emission asso- 
ciated with the shock surrounding the SW lobe, first discussed by 
Kraft et al. (2003), is X-ray synchrotron emission associated with 
particle acceleration at the shock front, and not thermal emission 
from shocked gas. We have investigated the spectral structure of the 
X-ray shell, presenting maps of spectral structure for the first time, 



and we conclude that non-thermal emission dominates over most 
of the shell, except in the inner regions, where the thermal X-ray 
emission can be used to determine the shock conditions. We have 
described a self-consistent model for the lobe and shock dynam- 
ics, considering in detail the implications of possible electron-ion 
non-equilibrium in the shocked material. Our analysis of the spa- 
tial distribution of power-law spectral index constrains ^y ma x for the 
accelerated particles to be ~10 8 at the strongest part of the shock, 
consistent with the expectations of diffusive shock acceleration the- 
ory. We considered the implications of our results for Cen A's role 
as a source of cosmic rays, but conclude that the inner lobe shocks 
are unlikely to be an important source of ultra-high energy cosmic 
rays as this would require very high magnetic field strengths. Fi- 
nally we considered the expected inverse-Compton emission from 
this synchrotron shell, finding that the predicted emission from IC 
scattering of starlight and the CMB is comparable to current upper 
limits at TeV energies. 
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